Rare Earth Monopnictides and Monochalcogenides from First 
Principles: Towards an Electronic Phase Diagram of Strongly 

Correlated Materials 

L. Petit/'^'Q R. Tyer,2 Z. Szotek,^ W. M. Temmerman,^ and A. Svane^ 

^Department of Physics and Astronomy, 
Aarhus University, DK-8000 Aarhus C, Denmark 
^ Daresbury Laboratory, Daresbury, Warrington WA4 4^D, UK 

Abstract 

We present results of an ab-initio study of the electronic structure of 140 rare earth compounds. 
Specifically we predict an electronic phase diagram of the entire range of rare earth monopnictides 
and monochalcogenides, composed of metallic, semiconducting and heavy fermion-like regions, and 
exhibiting valency transitions brought about by a complex interplay between ligand chemistry and 
lanthanide contraction. The calculations exploit the combined effect of a first-principles method- 
ology, which can adequately describe the dual character of electrons, itinerant vs. localized, and 
high throughput computing made possible by the increasing available computational power. Our 
findings, including the predicted "intermediate valent" compounds SmO and TmSe, are in overall 
excellent agreement with the available experimental data. The accuracy of the approach, proven 
e.g. through the lattice parameters calculated to within ^1.5% of the experimental values, and its 
ability to describe localization phenomena in solids, makes it a competitive atomistic simulation 
approach in the search for and design of new materials with specific physical properties and possible 
technological applications. 
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I. INTRODUCTION 



Increased computer power, together with rehable electronic structure codes, have opened 
up the possibihty of exploring electronic properties of huge numbers of materials.- A system- 
atic first principles based theoretical study of entire families of compounds has the potential 
of exploring new trends in their physical properties or discovering materials with improved 
characteristics of possible interest for technological applications. 

The ab-initio determination of materials properties requires electronic structure calcu- 
lations that are accurate and parameter- free. With the advent of density functional the- 
ory (DFT),- quantum mechanical calculations of the electronic structure of solids were 
put on a rigorous theoretical foundation. Nevertheless, the lack of knowledge of the ex- 
act exchange-correlation energy functional constitutes a substantial obstacle regarding the 
predictive power of the calculations, as different mechanisms are at work depending on the 
relative importance of band formation energy on one hand and on-site electron-electron in- 
teraction and electron localization on the other. Standard electronic structure calculations, 
based on the local spin density approximation (LSDA)^i^ to density functional theory, work 
very well for conventional materials like simple metals and their alloys, characterized by 
delocalized electron states and large band formation energies. However, for systems where 
localized electron states occur, such as the 4/'s in rare earth compounds, the LSDA suffers 
from a sizeable self-interaction (SI) error,- brought about by the local approximation to the 
exchange-correlation energy functional. Consequently, different assumptions and/or param- 
eters need to be invoked to describe the electronic structure of strongly correlated electron 
materials.- 

The self-interaction corrected^ local-spin-density (SIC-LSD) energy functional is free of 
the SI error and fulfils the requirement of being parameter free and of sufficient accuracy 
to study trends in properties of entire families of materials containing both localized and 
itinerant electrons.-"— Here we apply the SIC-LSD to provide global understanding of the 
ground state electronic structure of rare earth monopnictides and monochalcogenides and, 
through that, also demonstrate the predictive power and accuracy of this approach. Specif- 
ically, the present study covers 140 distinct compounds, RX, with R referring to the rare 
earth atoms Ce, Pr, Nd, Pm, Sm, Eu, Gd, Tb, Dy, Ho, Er, Tm, Yb, Lu, and X standing for 
the pnictide atoms N, P, As, Sb, Bi and the chalcogenide atoms O, S, Se, Te, Po. Most of 
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these compounds are known to exist in nature, crystallizing predominantly in the rocksalt 
structure. Due to the intricate dual nature of the 4/-electrons, as a function of chemical 
environment these compounds display an extraordinarily wide range of electronic, magnetic, 
optical, and magneto-optical properties with, among other things, potential practical appli- 
cations in the field of spintronics.— A considerable number of theoretical studies on the 
RX compounds have resulted in greatly improved insight into their properties.—"— However 
those calculations have often been guided by experiment, relying on parameters to repro- 
duce the observed properties. The goal of the present systematic study is to establish a 
paramater-free unified theoretical picture of the RX electronic properties and possibly pro- 
vide directions for further experimental investigations. With future applications in mind, 
we show how large scale computational infrastructures can potentially be exploited in the 
search for the novel rare earth related materials that play an increasingly prominent role in 
innovative technologies and the quest for renewable energies.— 



II. THEORETICAL BACKGROUND 



In the LSD approximation the total energy functional is decomposed as 

occ. 

j^LSD ^ 1^^) + u[n] + Ve,t[n] + E^!^[n^, nj, (1) 

a 

where T is the kinetic energy operator, and represents the interaction with an external 

potential. The electron-electron interaction is given as the sum of the Hartree term U[n\ 
and the local exchange-correlation term, E^^^[n^.,nt^. 

In the Hartree-Fock approximation, the expectation value of the electron repulsion oper- 
ator in some /" multiplet state \1/"(LS'J) may be written asi^>i^ 

{^''iLSJ) I V — I ^"(L^J)) = -nin - + aiE^ + aa^^ + a^E^ , (2) 

where "^"-{LSJ) is in general given by a linear combination of n x n Slater determinants. 
Since V — commutes with total L, S, and J, the above matrix elements are equal for Mj 

— rij 

quantum numbers belonging to one {LSJ) family. E^ are the Racah parameters,— which 
are related to the Slater two-electron integrals (defined later). The coefficients at depend 
on the particular LS'J-configuration chosen. 
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The quantity E'^ is the Coulomb repulsion between two /-electrons, and the Coulomb 
energy term in the LSD total energy functional ([1]) corresponds to the ^E^n^ part of the 
first term in Eq. (Ej). The term, —\E^n, in Eq. ([2]), cancels an equivalent contribution in 
the Hartree energy term, representing the interaction of an electron with itself. A similar 
cancellation does not occur in the LSD functional, which therefore carries an unphysical self- 
interaction.- This interaction, which tends to be insignificant for extended band state, may 
lead to uncontrollable errors in the description of atomic-like localized states, for example 
the /-electrons in the rare earths. Following the suggestion by Perdew and Zunger-, in the 
SIC-LSD methodology, the LSD functional (1) is corrected for this spurious self-interaction 
by adding an explicit energy contribution for an electron to localize. The resulting, orbital 
dependent, SIC-LSD total energy functional has the fornt^"^ 



where 



j^sic-Lsu ^ j^Lsn ^ _ ^^^^^^ 

occ. occ. 

Ai?... = = J2 [UK] + E^r[n^]} , (4) 

a a 
occ. 

Eso = J](^„|e(r)rs1^/.„). (5) 



The self-interaction energy (4) consists of the self-Coulomb and self-exchange-correlation 
energies of the occupied orbitals ipa with charge density tIq, and spin density = 
For itinerant states, the self-interaction 5f vanishes identically, while for localized (atomic- 
like) states S^^'" may be appreciable. Thus, the self-interaction correction constitutes a 
negative energy contribution gained by an electron upon localization, which competes with 
the band formation energy gained by the electron if allowed to delocalize and hybridize with 
the available conduction states. Different localized/delocalized configurations are realized by 
assuming different numbers of localized states - here /-states on rare earth atom sites. The 
SIC-LSD approach is fully ab-initio, as both localized and delocalized states are expanded 
in the same set of basis functions, and are thus treated on an equal footing. If no localized 
states are assumed, E^^^~^^^ coincides with the conventional LSD functional, i.e., the 
Kohn-Sham minimum of the E^^^ functional is also a local minimum of E^^'-^~^^^ . 

The spin-orbit interaction (5) couples the band Hamiltonian for the spin-up and spin- 
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down channels, i.e. a double secular problem must be solved. The spin-orbit parameter, 

e(r) = -1^ 

' c^dr' 

in atomic Rydberg units, is calculated from the self-consistent potential V. 

Of the remaining terms in Eq. ([2]), the term accounts for Hund's first rule. The E'^ 
term only contributes to the level spacing of excited multiplets and will not be relevant in 
a functional for the ground state energy. The E^ term, describes the level spacing between 
multiplets in the maximum spin configuration, i.e. this term accounts for Hund's second rule, 
by which the multiplet of maximum total orbital momentum, L, has the lowest energy. With 
respect to the average /" (maximum spin) energy (the grand bary center), the maximum L 
multiplet is lowered by -jE^, where j = (0, 9, 21, 21, 9, 0, 0) for n = (1-7), and n = (8-14) 
respectively. It is an atomic effect which results in an increased stability at 1/4 and 3/4 
filling of the 4/-shell, and this second Hunds's rule effect is often referred to as the tetrad 
effect (TE).— Although, both LSD and SIC-LSD take into account Hund's first and third 
rules, with respectively the exchange interaction and spin-orbit coupling included in the 
total energy functional, Hund's second rule is not defined for the homogeneous electron gas, 
underlying LDA.^- To account for it in our total energy calculations, we add (a posteriori) 
to the SIC-LSD total energy functional the relevant correction 

AEt = ~jE\ (6) 

The E^ parameter is equivalent to Racah's 5-parameter,— and is given in terms of the 
reduced Slater integrals as: 

= i(5F2 + 6F4-91F6), (7) 

where Fk = F^/D^. Here are numerical constants (D^=225, 1089, 7361.64 for k=2, 4, 
6) and the Slater integral F^ is defined through 

Jo Jo v>) 

Here, is the f-partial wave as calculated in the self-consistent crystal potential, and r< 
(r>) denotes the smaller (larger) of the variables ri and r2. In reality, correlation effects 
in the solid environment tend to reduce the multiplet energy level splittings. Hence our 
calculated TE values are on average 15-20% larger than their experimental counterparts.— 
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III. CALCULATIONAL DETAILS 



Given the total energy functional e^ic-lsd ^ ^j^g computational procedure is as for the 
LSD case, i.e. minimization is accomplished by iteration until self-consistency. In the present 
work, the electron wavefunctions are expanded in the linear-muffin-tin-orbital (LMTO) basis 
functions.— The atomic spheres approximation (ASA) is used, whereby the crystal volume 
is divided into slightly overlapping atom-centered spheres of a total volume equal to the 
actual volume. A known shortcoming of the ASA is that different crystal structures have 
different degrees of overlap of the ASA spheres resulting in substantial relative errors in the 
evaluation of the total energy. While this inhibits the comparison of energies of different 
crystal structures, when comparing the energies of different localization/delocalization sce- 
narios within the same crystal structure the ASA error is of minor influence. To improve the 
packing of the structure empty spheres have been introduced on high symmetry interstitial 
sites. Two uncoupled energy panels have been considered when constructing the LMTO's 
to ensure an accurate description of semicore states. The LS'-coupling scheme is adopted for 
the localized /-states by starting the iterations with Wannier states of appropriate symme- 
try. During iteration to self-consistency the symmetry of the Wannier states may change, 
however grossly retaining their overall characteristics due to the fact that the energy scale of 
spin-orbit interaction is smaller than that of exchange, but larger than that of crystal field 
for the /-states. The rocksalt structure with the magnetic moments arranged ferromagnet- 
ically was assumed in all the calculations. To conduct the present large scale systematic 
study of rare earth compounds, in addition to a reliable theory and computer codes, we 
have made an extensive use of grid computing and automated data/metadata management 
tools..2^'24 

IV. TOTAL ENERGIES AND VALENCY TRANSITIONS 

In the SIC-LSD methodology, both the itinerant and localized limits are described by the 
same energy functional, which enables us to predict the ground state of correlated electron 
materials from total energy considerations.-i^ Hence, one can investigate localization phe- 
nomena in solids^ by realizing and studying different localization/delocalization scenarios, 
giving rise to different valency configurations, and through total energy minimization de- 
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FIG. 1: a) Total energy as a function of lattice parameter for the divalent, trivalent and tetravalent 
configurations of the Ho ion in HoSe. b) The energy differences between the divalent and trivalent 
configurations, AEjj-m, for the rare earth selenide series (RSe). Here the squares refer to the 
calculations without the tetrad effect, while the stars give the results with the tetrad effect included. 



termine the ground state electronic structure and valency configurations of the compound 
under consideration.— This nominal valency is defined as an integer number of electrons 
available for band formation, 

N,al = Z- N,ore ' NglC, (9) 

obtained by subtracting from the atomic number (Z) the sum of core (Ncore), and localized 
(Nsic) electrons. In this paper we will be using two interchangeable nomenclatures, /" and 
R*""*", to describe the configuration of the rare earth ion, implying n = Nsic and m = N^ai, 
respectively. The total number of /-electrons may be larger than n, because, in addition 
to the n localized /-states, the band states also contribute to the total /-count on a given 
ion. Note that our calculated valencies, refer to the number of rare earth electrons that 
contribute to bonding, and thus do not necessarily coincide with the nominal (ionic) valency 
of a compound. 

In our study of the rare earth monopnictides and monochalcogenides, for each given com- 
pound the self-consistent calculations involve the total energy minimization as a function of 
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lattice parameter with respect to a number of different possible localization/delocalization 
scenarios, in order to determine the global energy minimum and the corresponding equilib- 
rium volume. To understand what is involved in this process, we explain the minimization 
procedure on the example of HoSe depicted in Fig. [T^. What is seen in the figure are total 
energy curves, as a function of lattice parameter, for three energetically relevant valency 
configurations. The way they come about is based on electronic structure considerations. 
The Ho atom contributes 67 electrons to the total electron count of HoSe, 54 of which are 
considered to be part of the core (and semicore), and thus not contributing signifcantly to 
its chemistry. The remaining 13 electrons on the outer most shells, 11 of which are the 
/-electrons, are chemically active and constitute the valence electrons. The latter can either 
hybridize with the other band electrons or localize on site, playing no role in band forma- 
tion. Assuming all the 11 /-electrons to be localized, means that only two s-electrons of Ho 
contribute to band formation, a scenario which we refer to as divalent and which represents, 
with respect to /-electrons, the fully localized limit. Configurations with one or two delo- 
calized (bonding) /-electrons are respectively referred to as trivalent and tetravalent. The 
LDA scenario, which treats all the ligand and rare earth valence electrons, including the 11 
/-electrons, as band electrons, is found energetically unfavourable by ~ 1 Ry, as compared 
to any of the other scenarios, and is therefore not shown in Fig. [1^. Altogether at least some 
forty self-consistent calculations are required to establish the global energy minimum for a 
compound like HoSe which, as can be seen in Fig. [TK, is obtained in the trivalent scenario. 

From the calculated local energy minima we can evaluate the energy differences between 
respectively the divalent and trivalent configurations, AEjj-m, and the tetravalent and 
trivalent configurations, AEjv-m, both of which are a measure of stability of the ground 
state configuration. From Fig. [1^, we find for HoSe AEjv-iii= 90 mRy and AEu~iii= 20 
mRy, indicating a rather stable trivalent ground state, with no tendency towards delocal- 
ization of an additional /-electron, but with some, albeit marginal, tendency to localization 
of an additional electron. Thus whilst we predict a trivalent ground state at ambient condi- 
tions, with respect to possible phase transformations, for example under negative pressure, 
AEjj_jjj is the relevant energy scale to consider in HoSe. The calculated trivalent equilib- 
rium lattice parameter is equal to 5.55 A, differing from the divalent minimum by about 
5%. 

In Fig. [T)3, we present the calculated total energy differences, AEjj_jjj, for the entire 
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rare earth selenide series (RSe) from Pr onwards. The two curves shown refer respectively 
to the SIC-LSD calculations with (stars) and without (squares) the TE. The positive values 
indicate a trivalent ground state, whereas the negative values refer to the divalent config- 
uration being energetically more favourable. One can see that without the TE, AEii_jij 
decreases monotonously from trivalent PrSe to divalent EuSe, and again from trivalent GdSe 
to divalent YbSe. The increasing trend towards divalency is associated with the enhanced 
localized nature of the /-electrons and can be understood as a manifestation of the well 
known lanthanide contraction. A small AEjj^jji indicates that for the given compound no 
clear ground state emerges, as the divalent and trivalent configurations are close to ener- 
getically degenerate, due to the competition between band formation and self-interaction 
correction (localization) energies. 

The overall trend of decreasing energy differences is still recognizable in both halves of the 
selenide series when the TE is included. However, due to the abrupt decrease in AEfi_jjj at 
DySe a characteristic hump, also observed in experiments on rare earths^^i^, develops in the 
later part of the series, implying that this compound is considerably less trivalent than one 
would expect from the SIC-LSD result prior to taking Hund's second rule into account. Even 
more important are the changes in predicted ground state configurations for ErSe and TmSe 
(as well as PmSe), indicating that the tetrad effect can bring about a qualitative change in 
compounds where the original energy differences are small. Otherwise, inplementing Hund's 
second rule has only a quantitative effect on the AEjj^ju of those compounds that from 
the outset are either solidly trivalent or divalent. 

Taking the TE into account, the divalent-trivalent energy differences, AEjj^jjj, for all 
the RX compounds (except for the rare earth nitrides and the Ce compounds) are collected 
in Fig. [21 The positive and negative values again indicate respectively a trivalent and 
divalent ground state configurations. In Fig. [2^, for each of the ligand ions the AEji_jji 
data are plotted as clusters of bars, representing the rare earth ions ranging from Pr to 
Yb. We notice that the total energy trends, observed for the selenide series in Fig. [lb, are 
reproduced in all the ligands, but overall, in the chalcogenides the divalent configuration 
becomes increasingly more favourable as a function of ligand size, whilst in the pnictides the 
trivalent configuration remains the ground state, even for the "late" rare earth Eu and Yb 
(with the exception of EuBi). The trivalent to divalent localization transition thus mainly 
pertains to the rare earth chalcogenides. In Fig. [2)3, the same total energy data are depicted 
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FIG. 2: a) The energy difference AEji-iii as a function of ligand ion for all the rare earths, b) 
The energy difference AEjj-jjj as a function of rare earth ion for all the ligands. 



as a function of rare earth, with the hgands plotted as clusters of bars ranging from P to 
Po for each of the rare earth ions. In this representation one can clearly see the combined 
effect of the lanthanide contraction, leading gradually towards the more preferred divalent 
configurations, and the repetition of this trend due to the half-filled shell effect. 

Although AEjj^jjj is the relevant energy scale for the large majority of the RX, there 
are a number of exceptions where the delocalized tetravalent scenario is energetically more 
favourable than the localized divalent scenario even though according to our calculations it 
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FIG. 3: The valency phase diagram for the entire manifold of rare earth monopnictides and 
monochalcogenides. The blue areas mark divalent compounds, while trivalency is indicated by the 
yellowish colour. The red parts of the diagram refer to the close to tetravalent compounds (mostly 
Ce compounds). The light blue/white regions represents the border line cases between divalent 
and trivalent configurations. The colour scheme refers to the valencies, obtained by scaling the 
calculated energy differences such that a value of 3.5 indicates trivalent-tetravalent degeneracy as 
observed in CeN, whilst values smaller than 2.5 refer to divalent ground state configurations. 

never actually becomes the true ground state. The trivalent-tetravalent energy difference, 
AEjv-in, is the more relevant energy scale for the Ce compounds for example, as the 
divalent configuration is always very unfavourable, and in most cases could not even be 
converged. In the case of CeN, and to some degree in CeO, the tetravalent and trivalent 
scenarios are actually close to energetically degenerate. Since the divalent scenario is so 
highly unfavourable, the Ce pnictides and chalcogenides, as well as all the other rare earth 
nitrides, are not shown in Fig. [2J For the LuX compounds AEii_iii is not adequate, as a 
divalent scenario would here imply 15 localized /-electrons. 

As shown in Fig. [3l a unified picture for all the RX compounds emerges when we map the 
calculated energy differences into ground state valencies. Here the valency scale is defined 
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such that it runs from 3.5 for the tetravalent/trivalent degenerate compounds {Ejv — Eiii=0) 
to 2.5 for the trivalent /divalent degenerate {Ejj — Ejjj=0) compounds. The yellowish areas, 
that dominate throughout the phase diagram, demonstrate the overall preferred trivalent 
ground state configuration for the large number of the RX compounds, in particular the 
monopnictides. The trend towards delocalizing an additional electron is only noticeable 
for the CeX systems and some of the nitrides, as indicated by the more reddish areas of 
the phase diagram, whilst the trivalent to divalent localization transition occurs around the 
regions indicated in white. The blue zones represent the stable divalent configuration, which 
is the dominating ground state for many of the "heavy" chalcogenides, in combination with 
the rare earths situated in the middle and later parts of the series. 

V. ELECTRONIC STRUCTURE 

A. Density of states: chalcogen versus pnictogen chemistry 

The observed total energy trends can be understood in terms of variations in the ground 
state electronic structures brought about by changes in the atomic composition, the distance 
between atoms, and the relative position of the atomic energy levels. In the RX compounds, 
the electronic structure results mainly from the /- and d-orbitals on the R-sites overlapping 
with the p-orbitals on the ligand-sites, as well as with the /- and rf-orbitals on neighbouring 
R-sites. In figure HI using TbP (column 1) and TmTe (column 2) as representative examples 
for the pnictide and chalcogenide rare earth compounds, we depict the density of states 
(DOS) for three different valency scenarios, as realized respectively in the delocalized (Figs. 

and UK') , trivalent (Figs. andHb'), and divalent (Figs. Ht andHt') configurations. 

The LSD band picture, for both TbP (Fig. H^) and TmTe (jUi'), is characterized by a 
large, /-derived, DOS at the Fermi level. A noticeable difference is observed with respect 
to the relative strength of the exchange and spin-orbit interactions. In TbP the exchange 
interaction dominates, which results in the splitting of the /-states into majority spin states, 
situated at the bottom of the phosphorus p-band, and minority spin states, situated at the 
Fermi level. In TmTe the spin-orbit interaction has become more important, splitting the 
thulium /-states into j = 5/2 and j = 7/2 bands, with the exchange interaction being 
responsible for the additional splitting within each of the two groups of bands. 
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FIG. 4: Left column: Density of states of TbP respectively for the a) LSD, b) trivalent and c) 
divalent Tb ion configuration. Right column: Density of states of TmTe respectively for the a') 
LSD, b') trivalent, and c') divalent configuration of the Tm ion. Density of states: total (dotted 
black), ligand p-projected (solid blue), rare earth d-projected (solid green), rare earth /-projected 
(solid red). The Fermi level is indicated by the dashed vertical line. 



In spite of its interesting features, the band picture (Figs. UK and H^') is an unphysical 
representation of the RX compounds, as the large DOS at the Fermi level is not observed 
experimentally. The trivalent scenarios of Figs. |4)d and |4)d' are obtained by localizing 
respectively eight /-electrons in TbP (Tb'^+ = Tb(/®)), and twelve /-electrons in TmTe 
^'jj^3+ = Tm(/^^)). Compared to the LSD scenario the trivalent configuration is found 
to be energetically more favourable by ~ 700 mRy for TbP and ~ 1200 mRy for TmTe. 
In TbP the seven majority spin states , as well as one of the minority spin states, are no 
longer available for band formation and hybridization and have vanished from the band 
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picture in Fig. |3}o. Similarly in the DOS for TmTe (Fig. Hb'), the localization of the six 
j = 5/2 /-electrons as well as six of the j = 7/2 electrons, implies that only two /-states 
are available for band formation. Notice that in the DOS plots only the band states are 
displayed, i.e. the itinerant valence states, including the delocalized /-states. The localized 
/'s are not shown since the SIC-LSD approach, which after all is a one-electron ground state 
theory, does not give accurate removal energies of localized states due to electron-electron 
interaction effects,— and the neglect of screening and relaxation contributions.— 

The trivalent DOSs, in respectively Figs. andHb', differ noticeably from each other, 
as in the former the Fermi level is situated in the gap between the filled phosphorus p-band 
and the unoccupied /rf-band, whilst in the latter the Fermi energy is pinned to the sharp 
/-peak above the conduction band minimum. The P atom has three unoccupied p-orbitals 
and in the TbP crystal it can accommodate three electrons from the Tb ion through charge 
transfer and hybridization which results in the filled p-band of Fig. |4)d. The Te atom on the 
other hand has one additional p-electron in its outer shell, which is why the p-band in TmTe 
can only accommodate two electrons, and as a consequence of this, here in the trivalent Tm 
scenario, the third electron is forced to occupy the /(i- conduct ion band. 

The variation in p-orbital occupation is fundamental to the overall difference that we 
observe in figure [3] with respect to the ground state valency configurations of the pnictides 
and chalcogenides. Thus we find that when trying to localize an additional /-electron in 
TbP, the resulting divalent configuration, depicted in figure Ht, is energetically unfavourable 
compared to the trivalent scenario in figure Hb by ~ 250 mRy. The /® to /^ localization 
is associated with charge transfer, and brings about the depopulation of the P p-states, as 
indicated by the Fermi level being situated in the p-band. The removal of the bonding 
states results in a considerable loss in hybridization and Madelung energy, which is not 
compensated by the gain in self-interaction energy, and the divalent scenario will therefore 
not be realized in TbP. In TmTe on the other hand, the localization transition from /^^ in 
Fig. Hb' to /^^ in Fig. Ht', is associated with a minor loss in band formation energy, given 
the vanishing width of the occupied /-peak, and the gain in self-interaction (localization) 
energy is comparatively larger by ~ 35 mRy. Accordingly, the divalent scenario of Fig. Ht', 
becomes the ground state scenario for TmTe. 
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B. Electronegativity and lanthanide contraction 

Apart from the fundamental difference between the pnictide and chalcogenide electronic 
structure, the observed energy trends in each series, as a function of R-ion in figure [2^ 
and as a function of ligand in figure [2)3, are determined by further factors such as ligand 
electronegativity, and lanthanide contraction. These changes in electronic structure are 
presented schematically in Fig. |5l Here we depict the densities of states (DOS) for eight 
distinct scenarios observed throughout the rare earth compounds. Each DOS is composed 
of several rectangles, referring respectively to the ligand p-states (situated lowest in energy), 
the rare-earth (i-states, and the localized and delocalized /-states. The occupied states are 
indicated in blue. For both the pnictides and chalcogenides, the electronegativity decreases 
with increasing anion size, i.e. respectively from N to Bi and from O to Po. On the other 
hand, the /-electron localization increases due to the lanthanide contraction, as we move 
through the rare earth series from Ce to Eu, and again from Gd to Yb. In terms of DOS, 
these two trends manifest themselves respectively via ligand p-bands moving up in energy 
(towards lower binding energies), and unoccupied /-states moving towards higher binding 
energies. In Fig. \5[ the combined effect of these two trends is depicted in the two sets 
of DOSs respectively along the diagonals of the pnictide ((A) — t- (C)) and chalcogenide 
compounds ((D) — )■ (F)). As the /-states move towards lower energy they get pinned to 
the Fermi level and start filling up. Ultimately the degree of /-band filling becomes such 
that the gain in energy associated with localizing the given /-state becomes more important 
than the corresponding loss in band formation energy, and the divalent scenario becomes the 
ground state (scenarios (C) and (F')).— As seen in Fig. \5[ owing to the additional electron, 
in the chalcogenides the filling of /-states sets in much earlier than in the pnictides, which 
explains why the trivalent scenario is relatively more dominant in the pnictides. 

VI. ELECTRONIC PROPERTIES 
A. Lattice parameters 

For each given compound, the calculated global energy minimum determines the equi- 
librium lattice parameter as well as the ground state. The calculated lattice parameters 
for the rare earth sulphides and the ytterbium pnictides and chalcogenides are compared to 
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FIG. 5: Schematic representation of densities of states (DOS) for the rare earth monopnictides 
(upper panel) and monochalcogenides (lower panel). Each DOS is composed of a number of 
rectangles representing respectively the ligand p-states (situated lowest in energy), the rare earth 
(i-states (situated highest in energy), and the occupied and unoccupied /-states (upright bars). 
Occupied states are shown in blue. For each panel the three different scenarios ((A) to (C)) and 
((D) to (F)), represent the different trivalent DOS resulting from the decrease in electronegativity 
and the increase in /-state localization as one moves along the diagonal of the corresponding 
compound space. Scenarios (C) and (F') refer to the divalent scenarios obtained respectively form 
(C) and (F) through localizing an additional /-electron. 



their experimental counterparts in respectively Figs. ^ and Eb- The agreement is good, 
reproducing the sudden increase in lattice parameter at SmS, EuS and YbS, due to the 
increased /-electron localization. We have observed that quite systematically the calculated 
lattice parameters for the divalent and trivalent scenarios of a given compound differ by 
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FIG. 6: Lattice parameters of a) the lanthanide sulphides (RS) and b) the Yb pninctides and 
chalcogenides (YbX) as calculated in their ground state configuration. 

about 5 %. Therefore, a comparison to the experimental lattice parameters is expected to 
provide a good indication as to whether the calculations predict the correct ground state 
and valency configuration. In Fig. [TJ the experimental equilibrium lattice parameters for all 
the rare earth monopnictides and monochalcogenidea^ are plotted along the x-axis, whilst 
our calculated equilibrium counterparts are plotted along the y-axis. In this representation, 
the proximity of the data points to the x = y line is a measure of the agreement between 
theory and experiment. We find that for the large majority of the compounds the calculated 
equilibrium lattice parameters lie within ~1.5% of the experimental values, i.e. well below 
the 5% that would indicate a wrongly predicted ground state configuration. The only ex- 
ceptions are CeN, DyTe, and three of the polonide compounds. With respect to the latter 
compounds, polonium is difficult to handle, and it is unclear to what degree the experimental 
results refer to stoichiometric samples.— 

It is well established from rare earth halide datai^ as well as XPS measurements on the 
rare earth elements^^ that the divalent configuration becomes more competitive in Dy com- 
pounds due to the 3/4 filling of the /-shell. The inclusion of the TE leads to a considerable 
reduction of Ejj — Em for DyX compounds, as can be seen on the example of DySe in 
Fig. [T|d, and quite generally gives rise to the negative dip around Dy in Fig. [2]3. This 
trend reproduces overall very well the experimentally observed data,— however we wrongly 
predict a Dy^~^ ground state in DyTe. This seems to indicate that in some cases we tend to 
overestimate the tetrad effect's influence as a result of approximating the respective Slater 
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integrals by the fully localized atomic limit instead of considering also corrections due to 
the solid state environment (screening/hybridization). For CeN, the experimental lattice 
parameter seems to suggest a tetravalent ground state,— whereas our calculations predict 
the trivalent configuration to be marginally more favourable than the tetravalent configu- 
ration by ~3 mRy, which on the other hand could be indicative of an intermediate valency 
scenario implied by experiments.— 

In Fig. [3, the wrongly predicted ground state of DyTe shows up as noticeably off diagonal 
data point, but at the same time it also highlights the excellent agreement obtained for all 
the other cases investigated with the present approach. Thus, the SIC-LSD methodology is 
capable of predicting, with considerable accuracy, the ground state valency and lattice pa- 
rameter of the rare earth monochalcogenides and monopnictides, and the associated ground 
state electronic structures can be expected to make valid predictions concerning the other 
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FIG. 8: A global picture of the electronic properties of all the studied pnictides and chalcogenides, 
based on the densities of states at the Fermi energy Fjp. The colour scheme is as follows: metals are 
marked in blue tones, with light blue for semimetals. Semiconductors are in black, whilst yellow to 
red is assigned to increasingly heavy fermion-like systems. The fields marked with "IV" indicate 
intermediate behaviour between two different valency configurations. 

physical properties of these materials. 

B. Electronic Phase Diagram 

Here we concentrate on establishing global trends in the ground state electronic structure 
which we base on the calculated DOS at the Fermi energy, Ep, as providing a guideline for 
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predicting materials with specific properties when we systematically scan the full manifold 
of rare earth compounds. The corresponding results are summarized in the form of an "elec- 
tronic phase diagram" presented in Fig. [HI where based on the calculated DOS we distinguish 
between semiconductors, semi-metals, metals, and heavy fermion-like compounds. 

Looking at the phase diagram, we see that our calculations predict a considerable num- 
ber of compounds to be semiconductors. In particular, with respect to the pnictides, NdN, 
PmN and the majority of the rare earth nitrides from GdN onwards are found to be small 
gap semiconductors, with energy band gaps ranging from 0.06 eV in NdN to 0.6 eV in 
HoN. Few conclusive conductivity measurements for the nitrides exist, and it is for ex- 
ample still not fully established whether GdN is semiconducting^^ or semi-metallic— Our 
calculations predict GdN to be a semiconductor, although in an earlier SIC-LSD calcula- 
tion, where spin-orbit coupling was not taken into account, a half-metallic ground state was 
observed instead.— All the divalent monochalcogenides turn out to be semiconductors, in 
good agreement with experiment. The divalent character implies that in the chalcogenide 
semiconductors the energy gap is brought about by localizing an additional /-electron, as 
depicted through the transition from (F) to (F') in Fig. [5l in variance to the energy gap 
observed in the trivalent nitride semiconductors, resulting from the filling of the unoccupied 
anion p-states, as represented by scenario (A) in Fig. [51 In the SIC-LSD the occupied 
4/-levels tend to be situated at too low energies,— as this ground state methodology can 
not reproduce the spectroscopy of an eventual experimentally observed mid-gap /-level. For 
EuO, EuS, and YbS, the calculations overestimate the energy gaps by a factor 2.—'^ 

The dominant part of the electronic phase diagram in Fig. [H] consists of metallic com- 
pounds. The indicated differences in conductivities can be straightforwardly related to the 
trends in electronic structure observed when moving along the diagonals in Fig [5] In the 
pnictides, the closing of the semiconducting gap due to the onset of pnictogen p - rare earth 
sd overlap gives at first rise to a vanishingly small DOS, resulting in the semi-metals CeX to 
PrX and GdX to DyX (light blue areas). As the overlap gradually increases (due to decreas- 
ing electronegativity and/or increasing lanthanide contraction), the behaviour becomes fully 
metallic (dark blue areas), and eventually in the "late" rare earths, with more and more of 
the /-states getting occupied, heavy fermion-like behaviour sets in (yellow-red areas). The 
rare earth chalcogenides display less variety in terms of metallic behaviour. Here we start 
out with a partially occupied conduction band, and semi-metallic behaviour is not observed. 
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Instead, as can be seen from Fig. [HI the rare earth systems, CeX to PmX, and GdX to 
ErX (apart from the polonides and a couple of tellurides) are metaUic. Compared to the 
pnictides, in the late rare earth chalcogenides, the /-states start filling even before they 
have moved far enough towards lower binding energies to achieve significant hybridization. 
Instead of forming a heavy metal, a relatively pure /-state pins the Fermi level (scenarios 
(E) and (F) in Fig. [5]) and the divalent semiconducting scenario (scenario (F') in Fig. [5]) 
correspondingly becomes energetically more favourable. 

The agreement between observed and predicted properties is overall very good, except 
maybe for the calculated heavy fermion pnictides. Thus the "late" rare earths, SmX, and 
YbX, are described as low carrier semi-metallic compounds, seemingly at odds with the 
large DOS at the Fermi level that is predicted in the SIC-LSD calculations. The experimen- 
tally observed large specific heat coefficient in these systems has been interpreted as heavy 
fermion behaviour by some authors,— whilst others have explained it in terms of the Kondo 
effect.—*^ Since our calculations are based on the effective one electron picture, they can not 
adequately describe the many-body physics of these compounds. However the large peak at 
the Fermi level that is calculated in the trivalent configuration indicates strong ^/-mixing, 
i.e. a prerequiste for the hybridization between the narrow /-peak and the sea of conduction 
electrons that underpins the Kondo effect. EuSb, YbBi, EuBi (as well as EuAs) do not ac- 
tually crystallize in the NaCl structure investigated here,— i.e., the very same pnictides that 
we find to be divalent or close to divalent do not occur in nature. As indicated schematically 
in Fig. [5l for the pnictides, localization of an additional /-electron ((C) — (C) transition) 
implies the removal of some bonding p-states. This, however, would destabilize the under- 
lying NaCl crystal structure and, as a result of this competition between band formation 
and correlations, it seems that these materials prefer to distort and crystallize in a different 
structure altogether. 

In Fig. [8]the compounds marked as "intermediate valent" (IV) are those where the calcu- 
lations find close to degenerate energy configurations, meaning either AEu^jjj or AEjj^m 
is smaller than 5 mRy. Given the dual nature of the /-electrons, it is to be expected that 
these materials are highly sensitive to environmental changes such as pressure^^^ or dop- 
ing. In the phase diagram of Fig. [3] they correspond to the trivalent-divalent degenerate 
compounds in and around the white area, and to CeN, the only compound that is at the 
borderline of the tetravalent-trivalent degeneracy. For the Tm chalcogenides, it emerges 
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from XPS measurements that TmS is metallic trivalent, TmSe is intermediate valent, and 
TmTe is semiconducting divalent.— This, as can be seen from Fig. |8l exactly matches 
the valency sequence that we have derived from the SIC-LSD calculations. For SmS we find 
AEjj_iii = -5.6 eV, indicating a divalent semiconductor, but given the small energy differ- 
ence, susceptible to undergo a localization-delocalization transition under pressure,— as has 
been observed experimentally.— Apart from EuO, the rare earth monoxides do not appear 
to occur naturally, but a number of them have been synthesized under high pressure, where 
it has emerged that EuO and YbO are divalent semiconductors, CeO, PrO and NdO are 
trivalent metals, whilst SmO is a metal with the Sm in an intermediate valent state.— These 
findings, including the intermediate valent nature of SmO, are again in excellent agreement 
with our predictions. 

In summary, this effective one-electron SIC-LSD methodology accurately predicts the 
phase transitions associated with the /-electron localization/delocalization. However, we 
should note here that because of the many-body effects, beyond those included in the 
SIC exchange-correlation energy functional, the representation of the calculated densities 
of states in terms of the described above phase diagram is most likely not the full picture of 
the underlying electronic structure. This is in particlular true in the TV and heavy fermion 
areas of the phase diagram, where those 'left out' many-body interactions will introduce ad- 
ditional refinements into the physical trends based on SIC-LSD. Therefore, compounds that 
are situated in those areas provide the greater challenge and should be further investigated 
by methodologies such as the dynamical mean field theory (DMFT)^ or GW.— 

VII. CONCLUSIONS AND OUTLOOK 

We have shown that the orbital dependent DFT approach, based on the LSD energy 
functional, appended by two important corrections, namely the self-interaction correction 
and the tetrad effect, can provide a global picture of the electronic structure of the rare 
earth monopnictides and monochalcogenides, expressed through the uncovered trends in 
their behaviour. We believe that this SIC-LSD based study provides an important under- 
standing of the physical properties of these compounds. Although for a given compound the 
DMFT+LDA or GW approaches can give added insight into the details of the electronic 
structure, they can not at the present time be employed for large scale predictive studies 



22 



given that they either depend on parameters and/or are very time consuming. 

Our study has given rise to a unified picture of the ground state valency and valency tran- 
sitions across the entire range of rare earth monopnictides and monochalcogenides, starting 
from tetravalency, through trivalency to divalency, as one moves from the early to late 
rare earths, in correlation with the ligands. We have also discovered the important links 
between the underlying electronic structure, valency and physical properties of these com- 
pounds, ranging from semiconducting to semi-metallic, metallic and even heavy-fermion-like 
behaviour. 

The present study has also established that if the energy differences between the ener- 
getically relevant valency scenarios are small, then fiuctuations will play an important role, 
giving rise to intermediate ground states. This effect can be described within the SIC-LSD 
formalism by mapping the valence (and also spin) fluctuations on to disorder in the spirit 
of the Hubbard III approach,— which we intend to apply in future to the "intermediate" 
compounds identifled in the present study. 
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